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The flavour and quark mass dependence of thermodynamic quantities in 
lattice QCD * 

A. Peikert with F. Karsch, E. Laermann 

Fakultat fiir Physik, Universitat Bielefeld, 33501 Bielefeld, Germany 

The dependence of thermodynamic properties of QCD on the number of quark flavours is investigated. Lattice 
results for the equation of state are presented for 2, 2+1 and 3 quark flavours. The simulations have been 
performed with the improved p4-staggered fermion and a Symanzik improved gluon action on lattices of size 
16 3 x 4 and 16 4 . 

1. INTRODUCTION 

In this paper thermodynamic observables like 
the critical temperature T c /^/a and the pressure 
p/T 4 have been calculated in full QCD with an 
improved staggered fermion action. The use of 
improved actions in finite temperature calcula- 
tions is strongly suggested by results from the 
pure gauge sector which have shown that im- 
proved actions reduce finite cut-off effects in ther- 
modynamic quantities significantly 0,0. In ad- 
dition, the calculation of the pressure with the 
standard staggered action did show strong cut-off 
effects when comparing results obtained on lat- 
tices with temporal extent N T = 4 and 6 H . 
The relevant number of flavours in finite temper- 
ature QCD is expected to be two light and one 
heavier quark. To quantify the effect of vary- 
ing the number of quark flavours simulations with 
Nf ~2, 2+1 and 3 are performed. 




6 8 10 12 14 16 

Figure 1. The fermion free energy in the ideal 
gas limit. 

Fat-links lead to an improved flavour symme- 
try [|| which results in a reduced pion splitting 
measured by the quantity 



2. THE SIMULATION 

The fermion fields have been simulated using the 
tree-level p4 action [Q including fat-links ||] with 
a fat-weight lj = 0.2 in the one-link derivative. 
This action is constructed to improve the rota- 
tional symmetry of the free quark propagator up 
to 0(p 4 ). The high temperature ideal gas limit 
of the free energy is approached much faster for 
the p4 than for the staggered action (Figure 1). 
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Iii Figure 2, S is shown for the staggered and p4 
fat action measured at different quark masses. 
Comparing results at the same lattice spacing 
a one finds a reduction of the pion splitting of 
about a factor of 2. 

For the gauge fields the tree-level 1x2 action 
has been used which in pure gauge simulations 
lead to an improved finite cut-off behaviour for 
quantities like latent heatjlj and pressure [§]. 
The gauge and fermion fields have been updated 
with the standard Hybrid R algorithm with a 
step size AT<m q a/2 and a trajectory length 
of 0.8. The lattice sizes are 16 3 x 4 and 16 4 . 
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Figure 2. The pion splitting S for standard 
(Nf = 2) and p4 fat improved (Nf = 3) actions. 

3. THE CRITICAL TEMPERATURE 
AND THE TEMPERATURE SCALE 

The critical temperature has been calculated for 
various quark masses for 3 flavour QCD and in ad- 
dition also for 2 and 2+1 flavours at a quark mass 
ofm g = 0.10. The pseudo-critical coupling j3 c was 
determined by the peak position of the suscep- 
tibility of the Polyakov-loop and the chiral con- 
densate, respectively. At the pseudo-critical cou- 
plings zero-temperature calculations on 16 lat- 
tices have been performed measuring the string- 
tension yfa and the meson masses mpg and my. 
The results for T c /-Ja axe plotted in Figure 3. 
The qualitative behaviour of the 2 and 3 flavour 
values for the standard and p4 action is quite 
similar; in both cases T c drops quite fast already 
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Figure 3. The critical temperature T c /*Jcr for 2, 
2+1 and 3 flavours for the p4 action. For com- 
parison also 2 flavour results obtained with the 
standard action Q are plotted. 



at large quark masses. For the improved p4 ac- 
tion T c /y/a for different number of flavours does 
only show a small difference. 
To set the temperature scale for the pressure 
calculations zero temperature simulations have 
been performed at a variety of (3 values for 2, 
2+1 and 3 flavours. The resulting string tension 
data have been fitted to an ansatz proposed by 
Allton (see Figure 4) 

(oV5)(/3) - R(J3)(1 + c 2 a 2 (/3) + c 4 a 4 (/3) + ...)/<*> 

with a? {13) = R(f3)/R(/3 ), the two- loop beta 
function R(f3) and the normalization /3q- 
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Figure 4. The string tension data and the fit to 
the data for the Nf = 2 + 1 case. 

4. THE EQUATION OF STATE 

The pressure p/T 4 has been calculated for 2, 2+1 
and 3 quark flavours. The masses of the light 
quarks are m/T = 0.4, the mass of the heavy 
quark is m/T = 1.0. From the difference of the 
gluonic part of the action on finite temperature 



and zero temperature lattices (S 1 
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the pressure can be calculated according to 
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In Figure 5 the result for the differences of the 
gluonic action densities are shown. One observes 
an increase for the maximum of that quantity 
corresponding to the increase of the number of 
flavours. The integration over these curves gives 
the pressure p/T 4 plotted in Fi gure 6. For larger 
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Figure 5. The action differences (S 1x2 )q — 
(S lx2 ) T for N f = 2, 2+1 and 3 for the p4 
fat improved action. 

temperatures the 2+1 flavour pressure is closer 
to the 3 flavour than to the 2 flavour pressure 
as expected from the high temperature Stefan 
Boltzmann limit (see arrows in Figure 6). If 
one normalizes the pressure to the appropriate 
ideal gas value one finds the same temperature 
dependence for 2 and 3 flavours for the p4 ac- 
tion (Figure 7). Figure 7 also shows that the 
systematic deviations from the Stefan Boltzmann 
limit obtained with different fermion actions are 
qualitatively in agreement with what has been 
calculated as the cut-off effect of the free energy 
at N T = 4 for the different actions (see Figure 1). 
The cut-off effects of the 2 flavour pressure for 
the p4-fat action and the standard staggered 
action are compared in Figure 8. As the pres- 
sure calculated with the p4 action seems to scale 
with the relevant number of degrees of free- 
dom we may expect to obtain an estimate for the 
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Figure 6. The pressure for Nt = 2, 
2+1 and 3 for the p4 fat improved action. 
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Figure 7. The pressure normalized to the ideal 
gas value p/psb for N T = 4. 
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Figure 8. The 2 flavour pressure for the standard 
action ( N T = 4, 6), the p4 fat action (N T = 4) 
and the rescaled continuum pure gauge result. 

remaining cut-off dependence by comparing with 
the appropriately rescaled continuum extrapola- 
tion for the pure gauge pressure. This is shown 
as the black curve in Figure 8. If this curve ap- 
proximates the continuum extrapolated pressure 
of full QCD at high temperatures the p4 action 
indeed reduces strongly the cut-off effects. A 
direct verification will require a simulation on a 
N T = 6 lattice. 

REFERENCES 



B. Beinlich et al., Phys. Lett. B390 (1997) 268 

B. Beinlich et al., Eur. Phys. J. C6 (1999) 133 

C. Bernard et al. (MILC),Phys. Rev D 
55(1997) 6861 



U.M. Heller et al., [hep-lat/9901010 
T. Blum et al., Phys. Rev. D55 (1997) 1133 
M. Liitgemeier, Ph.D. thesis, Bielefeld, 1998 



C. Allton , |hep-lat/9610016| 



